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CN ' Abstract 

o ' 

^vq , The advent of high-energy hadron colliders necessitates efficient and accurate compu- 

tation of multi-jet production processes, both as QCD processes in their own right and as 
backgrounds for other physics. The algorithm that performs these tasks and a brief numerical 
Q H ' study of multi-jet processes are presented. 

J-h ; 1 Introduction 

With the recommissioning of the Tevatron, and the forseeable commencement of physics at 
the LHC, the need for fast and accurate QCD calculations is now larger than ever. In this 
paper, we describe our efforts to arrive at such calculations. To set the stage: we have in mind 
the computation of QCD cross sections with many observed jets that, as usual, are modelled 
by assuming that each jet comes from a single fragmenting parton. With the large amount of 
energy available, the number n of jets can easily be as large as 8, thus requiring the computation 
of amplitudes with 10 or even more external legs. 

Although in principle straightforward enough, the usual techniques of evaluating Feynman 
diagrams and integrating the resulting cross section by Monte Carlo are in practice hampered 
by the computational complexity of the problem. The following obstacles can be recognized: 

1. The flavours of the initial partons are never detected, and in most cases (barring, say, b 
tagging) neither are the flavours of the final-state partons. In what follows we take all 
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quarks (u, d, s, c, and b) to be essentially massless, although a version of the code with 
massive fermions exists. In the case of flavour integration, to be discussed later on, we 
always treat quarks as massless. This means that in any given jet configuration there are 
usually very many contributing processes. Even enumerating these is a nontrivial task, 



and calculating each individual process cross section is even more so [16 



2. In addition to flavour, neither the colour nor the spin of any parton is observed. For 
an amplitude with p quark partons and q gluons this implies that in principle (6) p (16) 9 
contributions have to be added. It is true that, in particular, many colour configurations 
lead to zero amplitudes, but figuring out precisely which these are is very hard. 

3. Each individual amplitude, with specified flavours, colours and spins, contains very many 
Feynman diagrams. In Appendix B we give a recipe for determining the precise number 
of graphs, at the tree level |L7|]. Typical results are that the process gg — ► 8g is described 
by 10,525,900 diagrams, and the process gg — ► 2g3u3u by 946,050 diagrams. Inclusion of 
loop corrections worsens this dramatically, of course. 

4. Each amplitude peaks in complicated ways inside the momentum phase space. Straight- 
forward integration is therefore impractical, and one has to search for efficient mappings 
to do importance sampling in a multi-particle phase space. 

All these difficulties are addressed in this paper, and here we describe our solutions, in reverse 
order. 

4. The peaking structure of the amplitude is dealt with by our phase-space generating al- 
gorithm SARGE ||, 11 1 . This algorithm is tailored to the generation of so-called antenna 



structures. Let p\, . . . ,p n be the momenta of the n partons involved. The peaking be- 
haviour of the cross section for the purely gluonic process involving n gluons, gg — ► (n — 2)g 
is then dominated by the following antenna structure: 

l 

(P1P2)(P2P3){P3P4) ■ ■ ■ {Pn-lPn)(PnPl) 

and any permutation of labels. Since processes involving quarks do not show other domi- 
nant peaking behaviour than the above ones, we can cover the n-particle momentum phase 
space with good efficiency. More details about SARGE can be found in p2l 



3. Over the last years new algorithms, along with their implementations, for computing the 
tree-order scattering amplitude have been proposed ||, |l^, |j. These do not involve cal- 
culation of individual diagrams, but rather reorganize, in a systematic way, the various 
off-shell subamplitudes in such a way that as little of the computation as possible is re- 
peated. The improvement in computational efficiency of these algorithms is nevertheless 
dramatic, down from about n! to something like 3 n . In the algorithm suggested originally 
in reference §§ the scattering amplitude was computed through a set of recursive equa- 
tions derived from the effective action as a function of the classical fields. These classical 
equations represent nothing but the tree-order Schwinger-Dyson equations, a fact that 
was already emphasized in subsequent approaches GJ and will be illustrated below for the 
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special case of QCD. In fact, recursive techniques have already been used in the past to 
compute multi-gluon amplitudes p). 

2. The usual spin and colour summation is replaced by a Monte Carlo integration. In the 
purely gluonic case this has been shown to work rather well. It must be realized that 
since we replace the more usual sum over discrete colour and spin states by a continuous 
(Monte Carlo) integration, the variance of the cross section can indeed be expected to be 
smaller. 

1. We tackle the flavour combinatorics in the same spirit: we assign to each quark and 
antiquark a flavour quantum number, that is, we extend the definition of the external 
legs by a direct product with a flavour vector in the space of the 5 available flavours. By 
taking these vectors randomly we can perform a sum over flavours by Monte Carlo as 
well, and the only computational difference between the case of one single flavour and that 
of / flavours is simply a factor /, whereas the discrete- flavour-sum combinatorics would 
lead to a much bigger loss in speed. The possibility of coherent superpositions of different 
flavours might seem awkward but is in fact quite natural since all quarks are treated as 
massless and therefore the distinction between flavours is to a large extent arbitrary. The 
only place where flavours are not treated on an equal footing is in the structure functions 
that describe the difference in probability of picking out different flavours: by a judicious 
weighting of the flavour vectors for incoming quarks and antiquarks we can handle this as 
well. 

Considering the Monte Carlo treatment of spin, colour and flavour, the essential point is to 
realize that we are entitled to use any representation for the corresponding information, as long 
as we end up with the correct sums on the average. 

Before finishing this introduction we want to mention that, so far, we can calculate par- 
tonic cross sections. These are typically good for global quantities like total cross sections, pr 
distributions and the like. A fuller treatment involves coupling the generated events to fragmen- 
tation programs like HERWIG [1C]. These programs typically require additional information on 
the 'spanning of the colour string'. In Q a way to do this for processes with zero or one quark 



lines is described; in [18] we shall indicate how the present program can be adapted in a similar 
manner. 



2 The algorithm for jet production 

2.1 The scattering amplitude 

Our starting point are the Dyson-Schwinger (DS) equations, which give recursively the re-point 
Green's functions. These equations hold all the information for the fields and their interactions 
for any number of external legs and to all orders in perturbation theory. We restrict ourselves 
to tree level calculations, so we discard ghost fields. Also we consider all quarks as massless. 
The couplings between fields, relevant for QCD, are shown below. 
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Coupling of the fields in QCD 

The recursive content of the DS equation, for the gluon field for example, can be understood 
diagrammatically as follows: 



p p 

%P2 



Opi 




J3 pi 

# P2 



Recursion equation for gluons 

The figure shows that a subamplitude with an off-shell gluon of momentum P, has contribu- 
tions from three- and four-vertices plus a fermion-antifermion vertex. The shaded blobs denote 
subamplitudes with the same structure. So, in effect, this is a recursive equation which we can 
write down immediately as (suppressing the colour): 



A»{P) =Y i 8 P=Pi A»{p i ) + (i 9s ) n^(p 1 )7V(P2)<7(pi,p 2 ) 

1=1 P= Pl +p 2 



+ 



V^ x (P, Pl , P2 )A„(p 1 )A x (p 2 ) ( T( Pl ,p 2 ) 



P=Pl+P2 



M) J- G^A u ( Pl )A x (p 2 )A p (p 3 )a( Pl ,p 2 ) (1) 



-P=Pl+P2+P3 
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where V puX (P,pi,p2) and G^ uXp are the 3 and 4-point vertices and 



" p2 



is the gluon propagator. The symbol a(pi,P2) - we will call it the sign function - takes into 
account the Fermi sign and has the values ±1. More on the sign function will be presented later 
on. For a quark of momentum P (suppressing the colour again) 



(2) 



i=i 



where S is the propagator 
and for an antiquark 



P=P1+P2 
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Recursion equation for quarks 



b ~p2 



1=1 P=Pl+P2 



(3) 




where 



Recursion equation for antiquarks 



In order to reduce computational complexity, as we will discuss later on, we replace the four 
gluon vertex with a three- vertex by introducing an auxiliary field H^ u . We rewrite the part of 
the QCD Lagrangian that describes the four-vertex in terms of the auxiliary field as follows: 



C H = -gf abc A*AlHi u ' c - H^H^ a 



(4) 



The recursion for the gluons changes slightly, namely only the part with the four-vertex. Addi- 
tionally, we also have an equation for the auxiliary field: 
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A^P) = ^2S P=Pi A^ Pi ) + (ig s ) n W"(Pi)T^(K)^(PbK) 

8=1 P=p 1 +p 2 

+ ^ E ^ A (^Pl ! P2)^(p 1 )^ A (p 2 )a(p 1 , P2 ) 

^=Pl+P2 

- ( 5s ) £ X^A u ( Pl )H Xp ( J)2 )a(p 1 ,p 2 ) 

P=p 1 +p 2 

(5) 

ff^(P) = -^ £ X^A x (p 1 )A p (p 2 )a(p 1 , P2 ) (6) 

-P=Pl+P2 

where X^ uXp is the new H-gluon-gluon vertex: 

These four equations, namely (||),(|6|),(|2|),(||), represent off-shell subamplitudes that are the build- 
ing blocks of any process. They are used iteratively, combining two (or three) momenta, at each 
step, to build a subamplitude. The iteration begins with the initial conditions for the external 
particles. In particular for a gluon we have: 

KiPi) = <?\(Pi)8*H, i = l,-n (7) 

where a is the colour, a = 1, . . . , 8, e^(p) denotes the polarization vector and i denotes one the 
external gluons. For the quarks and antiquarks the iteration starts with 



4>k(Pi) 



u{pi)5 kki if i incoming 

u{pi)5 kki if i outgoing 
u(pi)S kki if i incoming 

(8) 

u{pi)8 kki if i outgoing 
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where k is the colour, k = 1,2,3. The next step is combining two of the external momenta, in 
all possible ways, in accordance to the Feynman rules, to compute the next subamplitude. The 
iteration goes through in the same manner, using the equations (||), repeatedly, each 

time combining the new momenta obtained, with the remaining of the external ones. After n — 1 
steps there is only one momentum left to be combined, so the last step gives us the amplitude 
A(pi,p2, ■ ■ ■ ,Pn)- We refer to reference 0] for all further details of the algorithm. 



2.2 Colour and helicity 



In order to have an estimate of the production probability, one has to sum over all colour 
and helicity configurations. Summation over colours is a delicate subject. If one performs the 
summation in a straightforward way then one has to consider something like 8 nfl x 3 Hq x 3 nq 
configurations for the re-parton scattering, where n g ,n q ,nq is the number of gluons, quarks 
and antiquarks respectively. In this section we show how this summation can be replaced by 
integration, which is then suitable for Monte Carlo computation. As a first step a simplification 



of the colour structure is possible by defining the following object | 13fl 



Gab 



£ 

a=l 



A, B = 1,2,3 



(9) 



where G a is the gluon field and all other indices have been temporarily suppressed. The new 
objects are of course traceless, 3x3 matrices in colour space. The interesting property of 
this colour representation is that it leads to a "diagonalization" of the colour structure of the 
three-gluon vertex. More specifically the colour part of the three-gluon vertex is now given by 



/ 



abc + a + b + c 
l AB l CD l EF 



--^{5 AD $CF 5EB 



$AF 8cB Sed) 



(10) 




Colour flow in the gluon 3- vertex, as represented in Eq.iU 



This colour structure shows the colour flow in the real physical process, where gluons can be 
represented by quark-antiquark states in colour space and their self-interaction, as given by 
Eq.([lO|), reflects the fact that colour remains unchanged on an uninterrupted colour line. The 
recursion equations that include the gluon, like Eq.([|), now are modified according to Eq.(^), 
to reflect the new colour structure. The full content of the recursion equations, including the 
colour structure as just described, is listed in the Appendix A. 

This new, simplified colour structure of the vertices, allows us now to take the next step in 
making the computation of the colour part of an amplitude more efficient, by ridding ourselves of 
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the summations, mentioned in the beginning of this section, and replacing them by integration. 
To this end, we assign to a fermion a complex vector za, where the index runs from 1 to 3, 
representing its colour content. These vectors parametrize the 5-dimensional representation of 
SU(3) on the sphere and are subjected to the constraint: 

\Z1 | 2 + | Z 2 | 2 + \Z3 f= 1 

In this space, integration is defined through the proper definition for the invariant group measure, 
[dz\. 



\i=l / i=l 

We can use a simple polar coordinates parametrization to represent these complex vectors: 

Zl = e^ 1 cos 9 
Z2 = e 1 ^ 2 sin 9 cos £ 
z 3 = e 1 ^ 3 sin 9 sin ^ 
< fa < 2?r ; 0<6»<f ; < ^ < f (H) 

and in terms of these variables, the invariant measure becomes: 

—3^Y[J d(j>i\ J 2 d9 J 2 d( cos 9 sin 3 9 cos £ sin £ 

The colour structure as described by (|10D, shows that gluons can be interpreted as quark- 
antiquark pairs. So we construct the following vector, appropriate for describing the colour part 
of a gluon, which has this colour structure: 

3 

ri a (z) = V2iJ2 z*(t a ) ijZj a = 1, . . . , 8 (12) 

i,i=l 

where t a are the Gell-Mann matrices. This vector is real, because (rj a )* = rj a due to the 
hermiticity of the Gell-Mann matrices, and is normalized as follows: 



/ 



[dz] r] a (z)r] b {z) = 5' 



ub 



As far as our recursive equations are concerned their structure remains unaffected and the 
only thing to be changed are the initial conditions : 

GUPi) = E G a (Pi) V a (z) = V6 (z iA z* B - l -5 AB ^ e»(P t ) 

TpA(Pi) = V 7 ^ u(Pi) Z iA 

ij A (Pi) = V3 u(Pi) z* A (13) 
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where as usual i = 1, . . . , n, A is the helicity and z± are the new continuous colour coordinates 
of the i-th parton. The constants in the front are normalizations. 

In the same spirit summation over helicity configurations of the external partons can be 
replaced by an integration over a phase variable. For example, for a gluon this is achieved by 
introducing the polarization vector 

<(p)=e l V(p,+)+e- l V(p,-) 

where <p is a random number. Then by integrating over </> we obtain the sum over helicities, 

A=± 

The same thing can be used for the helicity of quark or an anti-quark. For example, for the 
quark we have 

U(f> {p) = e^u + (p) + e-^u.(p) 
and when integrated over 4> it gives the sum over polarizations 

u x (p)u\{p) =i 

and the same for an antiquark. 

2.3 The Fermi sign function 

Since we are dealing with fermions, we must find a way to incorporate a sign change when we 
interchange two identical fermions in a process. To this end we use the binary representation of 
the momentum labels of the external particles (e.g. Pi — > (0001), P 2 — > (0010), P3 = Pi + P 2 — > 
(0011), etc.) So the sign relative to the permutation of two momenta, cr(Pj, Pj) is computed as 
an operation on the two binary strings representing those momenta, a(mi,rrij). The function 
that performs that operation is defined as: 

a( mi ,m 2 ) = (-l)^ mi ' m ^ (14) 

with 

X(m 1 ,m 2 ) = ^2rh u { ^m 2j - J 

i=n \j=l / 

A hat over the binary string means that this particular bit is 
external particle is a boson. 

2.4 Flavor treatment 

The classification of processes contributing to the production of n jets is simplified when certain 
symmetries are taken into account. For instance processes like gg — > uuuu and gg — > dddd 
can be grouped together in the general class with representative gg — > qqqq; obviously an extra 




(15) 

set to if the corresponding 
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factor /, representing the number of light flavours, has to be taken into account. We may call 
the representative 'a distinct processes' and attempt a complete classification of such processes. 
Is clear that the number of distinct processes for a given number n of jets will be a small fraction 
of the total number of processes; this has a major implication on the computational complexity 
of the problem, since we may have the total contribution by just computing the contribution 
of a relatively small number of distinct processes. In the following table we give the number of 
distinct processes, as defined in Appendix C, for / = 5 final state flavours, and / = 4 initial 
state flavours, where we used these numbers in order to be able to compare with ref. |IJ. 



# of jets 


2 


3 


4 


5 


6 


7 


8 


9 


10 


# of dist. processes 


10 


14 


28 


36 


64 


78 


130 


154 


241 


total # of processes 


126 


206 


621 


861 


1862 


2326 


4342 


5142 


8641 



Another way to tackle all these different processes is to introduce quarks (antiquarks) that 
are a mixture of different flavours. Thus, keeping with the spirit of colour and helicity treatment, 
we integrate over flavours instead of summing over them. This is done by attaching a vector / 
to each spinor describing a fermion, the components of which are random numbers between 
and 1: 

4> P = u(p) x / 

with 

/ Si \ 

_ h 

f = ^/N] h 

V fN, J 

where fi, f%, . . . , f^ f are random numbers and Nf is the number of flavours. We choose the 
vectors / in such a way so that they are normalized as follows: 

(fifj) = Sij 

where i,j = 1, . . . ,Nf. As far as the final state is concerned, all flavours are equally treated 
(massless quarks). In the initial state however, due to different structure functions, special care 
should be taken of different flavours by a suitable weighting of the initial conditions with the 
structure function appropriate for each flavour. 

In treating flavour as described so far a new 'definition' of the concept of distinct process 
emerges. Any process is now composed by three primary objects, namely g, q and q, where q(q) 
represents a coherent superposition of all flavour states. The possible initial states are nine, 
namely gg, gq, gq, qg, qg, qq, qq, qq, qq and a final state will be determined by fixing the number 
m of qq pairs. Obviously m should satisfy, 2m < n — q, with q counts the net quark (antiquark) 
content of the initial state, namely 



c n 



1 c, 



qq 



^99 ~ — ^qq — ^ ^gq — ^qg — ^qg — ^gq 

So the number of 'distinct processes' is now given by 

9k + 3 if n = 2k and 9k + 7 if n = 2k + 1 



qq 



(16) 
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which represents a further reduction compared with the number of distinct processes indicated 
in the previous table. 



3 Jet production rates 

The features discussed in the previous sections have been implemented in a FORTRAN code with 
which we are able to perform perturbative calculations in QCD. We can divide the structure 
of the calculation leading to jet production rates in four parts: event generation, phase-space 
generation, squared amplitude computation and sum over processes calculation. In somewhat 
more detail these parts consist of: 



Amplitude computation The basic structure of the algorithm for the amplitude was dis- 
cussed in section 2. In order to obtain squared amplitudes we use projection to a continuous 
colour space and colour MC integration and integration over helicity states. What remains to 
be done is the ability to extract information about the colour string of the final states so that it 
can be used in codes like HERWIG [10]. Work is being done towards this end [18]. 



Phase space generation In the current status of the code we have implemented two ways 
of phase-space generation. First of all we have RAMBO ||, a flat phase space generator which 
has been around for many years and has proved reliable for the relevant calculations. Recently 
though, a new generator has been developed, SARGE [11], which take into account the so called 
antenna structure of QCD amplitudes. So far it has been proved very efficient for the generation 
of this type of phase-space S. 



Event generation Let us consider the scattering of two hadrons. The content of a hadron is 
characterized by the parton structure function f(x,Q 2 ) where x is the fraction of the momentum 
P of the hadron, carried by the parton, p = xP, and Q 2 is the QCD scale. The cross section 
for the scattering of two hadrons is given by the sum of all subprocesses between the parton 
constituents of the hadrons, weighted with the corresponding structure functions of the incoming 
partons: 




where s is the centre of mass(CM) energy squared at the hadron level. As usual the functions 
F m (x) are defined as F m (x) = f m (x)/x, where f m (x) are the various parton structure functions. 
The is the matrix element squared, summed (i.e. integrated) over helicity and colour 

degrees of freedom, and d<& is the element of the phase space. The sum (i.e. integral) is over all 
partonic processes. We use Eq.(fl7|) to estimate jet production rates. 



Sum over processes Finally, once given the number of jets, we are able to compute all the 
relevant subprocesses that contribute at one go. This is done by randomly choosing a subprocess 
and then using Monte-Carlo to obtain the total cross section from all contributions. The random 
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choice of a subprocess is based on the choice of a pair of integers (i, m), i selecting one out of the 



nine possible initial states, and m being the number of qq pairs in the final state, see Eq.(l£). 

There is also the option to use an optimization based on the fact that some processes over- 
whelm the total cross section over others (for example the purely gluonic process, gg — > ng has 
the largest cross section by an order of magnitude, compared to processes with different initial 
or final states). 

The algorithm has been used firstly to compute cross sections for 4 and 5 jet production. 
We have chosen a CM energy of y/s = 14 TeV and we apply the following cuts: 

p T i > 60 GeV, i5 > 30° \rn\ < 3 



where pt = \lp% + Py is the transverse momentum of a jet, 6ij = Z. (j>i,Pj), is the angle between 

jets and rj = — In tan (|), is the pseudo-rapidity. 

For convenience all results are obtained with a non-running strong coupling constant put 
equal to unity There are several parametrizations for the parton structure functions in the 
literature. The one that we will use is the MRST parametrization |15| , and the number of light 
flavours is taken to be / = 5. 

In Table 1 we have listed several subprocesses and their cross sections, relative for pp scatter- 
ing at LHC. The table is organized as follows: we are not referring to a particular quark flavour, 
so we use the letters q, r, s to denote quarks or antiquarks. Since we consider all quarks and 
antiquarks massless it makes no difference which particular flavour appears in the final state. We 
compare our results with those from a well established code, NJETS. We also show distributions 
of the maximum px for two processes, namely gg — > gggg and qq — > rfss. 

To show the potential of the flavour integration method we calculate the total cross section 
for a set of possible final states and compare with NJETS where as usual a sum over different 
flavour states has been used. The results are presented in Table 2. We have used the letter q to 
denote any flavour and all combinations have been taken into account in the calculation of the 
cross section. We have used the same cuts and CM energy as before. 

We have also obtained the total cross sections for the production of 3,4,5,6,7 and 8 jets. 
These are listed in Table 3. Also shown there is the contribution of the purely gluonic process to 
the total cross section. We see the fraction of this process compared to the total is diminishing 
with increasing number of jets, as was already commented in ifjlj . 

Finally we have plotted distributions of the total cross section for the production of 5,6 and 
7 jets. The quantities shown are: 

• Transverse momenta pt of the products 



Invariant masses My = \Jl{pi ■ Pj) 



4 Summary 

We have presented a procedure to calculate matrix elements and cross sections in QCD using 
an iterative algorithm based on the Schwinger-Dyson equation. Thus we free ourselves from the 
task of computing all Feynman graphs for a process, a task which can become impossible even 
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for a moderate number of particles involved. We have also managed to bypass the computa- 
tionally expensive procedure of summing over all possible colour and helicity configurations, by 
introducing continuous vectors to represent these otherwise discrete quantities and using Monte- 
Carlo integration instead of summation. At this stage our code can reliably compute scattering 
amplitudes and cross sections 

• for individual, single-flavour processes 

• for processes where a sum over all contributing flavours is needed both in the initial or the 
final states 

• for total jet cross sections, where we are interested in production rates with contribution 
from all possible subprocesses and all flavours. 

Our future interests involve a convolution of this code with fragmentation codes like HERWIG, so 
that one can perform realistic simulations of multi-jet processes. 
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y/s = 14 TeV 


Process 


OUR CODE(n 


b) 


NJETS(nb) 


55 5555 


2.681 




2.533 




99 -»• 5555 


0.0020 




0.0021 




95 9555 


1.131 




1.159 




55 9955 


0.106 




0.104 




99 9955 


0.062 




0.059 




gg -> rf^fir 


0.586 x 10 


-3 


0.558 x 10 


-3 


99 9955 


0.134 




0.126 




gf — > gfgg 


0.171 




0.161 




gr — > grgg 


0.210 




0.197 




95 9995 


0.035 




0.033 




gfir _, rf gfir 


0.032 




0.035 




55 9999 


0.524 x 10 


-3 


0.526 x 10 


-3 


55 — ¥ QQ r ^ 


1.059 x 10 


-3 


1.074 x 10 


-3 


99 qqqq 


0.807 x 10 


-3 


0.851 x 10 


-3 


99 — * gg?~f 


0.866 x 10 


-3 


0.786 x 10 


-3 


gg — > rfrf 


7.94 x 10" 


-6 


7.93 x 10" 


6 


gg — > rf ss 


1.54 x 10" 


-5 


1.49 x 10" 


5 


99 9999 


1.64 x 10" 


-3 


1.60 x 10- 


3 


99 — * gg?*f 


1.67 x 10- 


-3 


1.60 x 10 - 


3 


qr — > grgg 


2.78 x 10" 


-3 


2.91 x 10- 


3 


gr — > grss 


2.63 x 10" 


-3 


2.57 x 10" 


3 


gf — ► gfgg 


2.28 x 10" 


-3 


2.41 x 10- 


3 


gf — > qfrf 


2.20 x 10" 


-3 


2.07 x 10" 


3 


qf — > gf ss 


2.39 x 10" 


-3 


2.12 x 10- 


3 



y/s = 14 TeV 


Process 


OUR CODE(pb) 


NJETS(pb) 


55 99999 


159.49 


160.40 


qq 99999 


0.141 


0.136 


qg qgggg 


87.94 


84.37 


gg qqggg 


9.193 


9.218 


qq ->• qqggg 


4.37 


4.561 


qq -> rfggg 


0.034 


0.0348 


qq qqggg 


11.47 


12.35 


qf -> qfggg 


11.39 


11.98 


qr -> qrggg 


16.51 


17.29 


qg qqqgg 


3.858 


3.799 


qg -> rfggg 


3.738 


3.887 


55 99995 


0.104 


0.101 


gg -> ggrfg 


0.207 


0.206 


99 99995 


0.252 


0.259 


99 — qqrfg 


0.230 


0.254 


qq — ► rfrfg 


0.0020 


0.0020 


qq — > rfssg 


0.0038 


0.0038 


99 99995 


0.684 


0.691 


gg — ► qqrfg 


0.698 


0.659 


qr -> grggg 


0.998 


0.922 


gr — > qrssg 


0.862 


0.941 


qf -> gfggg 


0.684 


0.658 


gf — > qfrfg 


0.650 


0.678 


qf — ► qfssg 


0.634 


0.665 


qg 99999 


0.0352 


0.0334 


95 — * qqrfq 


0.0646 


0.0682 


qg — > rfrfq 


0.0328 


0.0328 


qg — > rfssq 


0.0650 


0.0668 



Table 1 Production rates for 4 and 5 jet production. All results have 
an estimated error of 4% . 
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Figure 1: Maximum pt distributions for gg — ► gggg (top plot) and qq — ► rrss (bottom plot) 



= 14 TeV 


Process 


OUR CODE (nb) 


NJETS (nb) 


99 3955 


0.242 


0.232 




0.015 


0.014 


99 9955 


0.311 


0.322 


gg -> g<M 


0.013 


0.014 


95 9995 


0.167 


0.174 



Table 2 Flavour integration with an estimated error of 2%. 



# of jets 


3 


4 


5 


6 


7 


8 


a(nb) 


91.41 


6.54 


0.458 


2.97 xl0~ 2 


2.21 xl0~ 3 


2.12 xlO" 4 


% Gluonic 


45.7 


39.2 


35.7 


35.1 


33.8 


26.6 



Table 3 Total cross sections for the production of up to 8 jets. 
The last row shows the percentage contribution of the purely 
gluonic process ( gg — ► ng ). The estimated error is 3% for 3 
and 4 jets, 4% for 5 and 6 jets and 6% for 7 and 8 jets. 
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Minimum and maximum p T 



5 JETS 



1500 

p,!-« (GeV) 



6 JETS 



500 1000 1500 2000 2500 3000 3500 4000 4500 5000 
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Distribution of trail verse momenta 
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Figure 2: pt distributions for 5,6 and 7 jet production 
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.3 Distribution of invariant mass 
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;ure 3: Invariant mass distributions for 5,6 and 7 jet production 
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Appendix A: The complete recursion relations 

Below we give the full recursion equations including the colour. 
• Gluon field A"(p) 

A ab(P) = ^2 E {MPi)l^A{P2)-\{Y,MPi)l^c{P2)\5 AB }a{ Pl ,p 2 ) 

P= P1 +p 2 V C J 

+ ^2 E ^r( P .Pl'P2){^ C (Pl)^CB(P2)-^c(P2)^B(Pl)}^(Pl,P 

P=pi+p 2 

+ ^ E ^{^(Pl)^(K)-^(»)^ B (Pl)}^(Pl,fl2) 

P=r>i -l-r><-> 



P=pi+P2 

Auxiliary field H /xl/ (p) 

HTb(P) = ^ E ^{^(Pi)^(») " ^(R2)^b(pi)> "(pi.m) 

P=Pl+P2 

Fermion field ipA(p) 

M P ) = jk E ^4ab(Pi) ^B(P2) ff(Pl,P2) 



P=Pl+P2 



• Anti- fermion field i/)a{p) 



V^(-P) = J| E ' t l J B(p2)4BA(Pl)f °(Pl,P2) 
P=Pl+P2 

Appendix B: Counting diagrams 

In this appendix we show how the number of tree-level graphs for a general QCD process can 
be determined. This is done in the same recursive manner as that in which the amplitudes 
are computed. Let there be k quark flavours. We shall compute the number of graphs that 
take a single specified parton into a number of specified partons, as follows. Let us denote 
by ai(no, nl, ni,ri2,n2, ■ ■ ■ , nk,nk) the number of graphs with a quark of type i coming in and 
ending up in uq gluons, n\ quarks of type 1, n\ antiquarks of type 1, and so on. We define the 
following generating function: 



lfj(z,Xl,Xl,X 2 ,X2, ■ ■ ■ ,Xk,Xk) = 

V* a i (n ,n 1 ,n 1 ,n 2 ,n2,...,n k ,n k ) 
^ n !ni!ni!n2! • • -n k \ 

n ,n 1: ...,n k >0 



~n „ni-ni ri2 . . . ^n k 
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The generating functions tp for incoming antiquarks, and cp for incoming gluons, are denned 
in an analogous manner. The Feynman rules of QCD then tell us that the various generating 
functions are related in the following manner: 

tpi = Xi + ipi4> , 
i>i = xi + Ppi(f> , 



11 _ 



6 

These can trivially be solved in such a way that everything is expressed in terms of (p: 

Ipi = Xi/(1 - (p) , Ipi = Xi/(l - (p) , 

and 



k 



X{Xi . 



The number of different quark flavours is seen not to be a complication here. By computer 
algebra, the function (p can easily be obtained to quite high order in z and the x's, and the 
number of diagrams read off from the corresponding coefficient. A similar discussion can be 
found in (7j, but we feel that the above procedure is simpler and moreover leads itself to estimates 
of the asymptotic number of graphs for very large multiplicities [17]. 



Appendix C: Counting distinct processes 

We start by tabulating the five types of initial states (not counting trivial charge conjugation), 
with the corresponding number of distinct processes, and the functions needed to compute their 
multiplicity factors. We denote by n the number of final state partons: 



initial-state type 


distinct processes 


multiplicity factor 


A (gg) 


Ci(n) 


X(n ,ni,...,n f ;f) 


B (qq) 


C 2 (n) 


X(n ,n 2 ,...,n f ;f - 1) 


C (gq and qg) 


C 2 (n-l) 


X(n ,n 2 ,...,n/;/- 1) 


D (qq) 


C 2 (n-2) 


x(n ,n 2 ,...,n/;/-l) 


E (qq' and qq') 


C 3 (n-2) 


X(n ,n 3 ,...,n f ;f -2) 



In order to clarify what we mean we consider the example of the type A initial state. Each 
distinct process is defined by an array (no,Wi, • • • , Uf). For instance, in the case of four-jet 
production we have 

(4,0,0,0,0,0) gg^gggg 
(2,1,0,0,0,0) gg ggqq 
(0,2,0,0,0,0) gg qqqq 
(0,1,1,0,0,0) gg qqrf 
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Therefore, in order to count the distinct processes we need the following three functions: 
d(n) = 0(^i > n 2 > ... > n/) 

no+2ni+...+2n^=n 

C 2 (n)= G(n 2 > n 3 > ... > n f ) 

no+2ni+...+2rif=n 

and 

C 3 (n) = Yl e(n 3 >n 4 >...>n / ) 

no+2ni+...+2ny=n 

Of course each distinct process, given by the array (no,ni, • • • > n /) has a multiplicity factor 
that is easily computed: 

X(n ,ni,. . . ,n/;/) = n/(n/ - l)...(n/ - j + l)/j! 

where j is defined as 

/ 

j=z if n n ^° 

i=i 



i = / - 1 if nw ° 



j = l if m / o 

j = otherwise 

Results for / = 5 final state flavours, and / = 4 initial state flavours, are shown in the table 
of section 2, in the discussion of flavour treatment, just to compare with ref. |TJ. 

Moreover it is very easy to produce the list of the distinct processes to be computed for each 
case as well as the multiplicity factors x- A code doing this is available. 

To study high-n behaviour we may use the generating function technique. Then we get 

oo / 

f iW = y> ( „ K= n 



n=0 



Fo(x) 



(l- x ) 11(1 -xV) 
Y C 2 (n)x n = ————— {l _ x 2 j) 

n=0 v /v ; j=l y ' 

oo , f-2 



f 3 ( X) = E g 3Nx"= (1 _ x)( i 1 _ x2)2 n ir ^ X 2 J) 



n=0 



As one can easily see, that the order of the pole at x = 1 is always given by / + 1. 

Our results may be compared directly to those obtained in ref. |IJ. In fact, using our method, 
we were able to detect several errors in Table 9.3 page 125, among which the fact that a processes 
is missing for the cases m = 6 and m = 7, namely qr — > qrrf(g). Correcting for this error we 
get full agreement. 
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